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Abstract. We present a method to study phase transitions in the large N limit of 
matrix models using matched solutions of Whitham hierarchies. The endpoints of the 
eigenvalue spectrum as functions of the temperature are characterized both as solutions 
of hodograph equations and as solutions of a system of ordinary differential equations. 
In particular we show that the free energy of the matrix model is the quasiclassical 
T-function of the associated hierarchy, and that critical processes in which the number 
of cuts changes in one unit are third-order phase transitions described by matched 
solutions of Whitham hierarchies. The method is illustrated with the Bleher-Eynard 
model for the merging of two cuts. We show that this model involves also a birth of a 
cut. 
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1. Introduction 

We consider the partition function of the unitary ensemble of random Hermitian 
matrices [U |2] 



N 



/ + 00 p + oo IN \ 

■■■ - T s n(^^ - ^^•)' • • ■ d^^' 
■°° \ ^ i=l / i<,j 

where T > plays the role of a temperature, and where the potential V{z) is a real 
polynomial of even degree 

2p 

V{z)=J2tnZ-^ (2) 
n=l 

with V{0) = and positive leading coefficient t2p > 0. We are interested in the critical 
processes which arise in the large limit as T varies. It is known [3l H] that for a given 
temperature the support J of the eigenvalue density p{x) is the union of a finite number 
of real intervals 

s 

J=\J{P2,-l,P2j), /^l</32<---</92.. (3) 

The density and its support are uniquely determined by the following conditions [1]: 

/oo 
log \x — y\p{y)dy = L, for x G J, (4) 
-oo 

/oo 
log \x — y\p{y)dy > L, for x ^ J, (5) 
-oo 

for some real constant L, together with the normalization 

J p{x)dx = 1. (6) 
These conditions mean that the infimum of the functional 

/oo /"OO /'OO 

V{x)(f){x)dx — T / / \og\x — y\(j){x)(f){y)dxdy, (7) 
-oo J — oo J — OO 

on the set of normalized functions 0(x) > in (— oo, oo) with bounded support is 
attained at = p (the constant L is the Lagrange multiplier). In physical terms p{x) 
represents the equilibrium density of a normalized one-dimensional charge distribution 
in the presence of an external electrostatic potential V{x). 

The general picture of the model is that the T > semi-axis decomposes into a 
number of disjoint (not necessarily connected) subsets Ig wherein the support of the 
eigenvalues has a fixed number s of disjoint intervals. The multicut case s > 1 is 
essentially different from the one-cut case: in the latter, the partition function has a 
regular large N topological expansion, whereas for s > 1 no such regular expansion exists 
and oscillatory terms appear. Using a saddle-point argument it has been shown [5] that 
the asymptotic behavior of the partition function for s > 1 is of the form 

Zjv ~ e~^'^e{N/2'Kn), (8) 
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where 6{z) is a Riemann theta function and 

F = —T / V{x)p{x)dx + / log \x — y\p{x)p{y)dxdy. (9) 

J J JjJj 

In this paper we propose a method to study phase transitions Is — )■ Is±i such 
as the merging of two cuts [H [71 [8] or the birth of a cut [9l HO]. Our method uses 
techniques of the theory of quasiclassical integrable systems [H] (universal Whitham 
hierarchies) to characterize the endpoints /3j of the support as functions of T and 
of the coefficients t„ of the potential V. Following Krichever's algebro-geometric 
scheme [TT|, [121 ^Si, we derive hodograph type equations for the functions (3j (which 
are determined from the expressions of the Abelian differentials of an underlying 
Riemann surface) and characterize the free energy ([9]) as a quasiclassical r-function 
Although our approach is much inspired by well-known works on the analysis of 
the zero-dispersion limit of the Korteweg de Vries equation in terms of Whitham 
equations [Tl|15l|16l[I71[l8l[l9l[20l[2ll[22], we provide an additional efficient tool 
to study solutions of Whitham equations in matrix models: we prove that the /3j satisfy 
a system of first-order ordinary differential equations which is completely determined by 
the expressions of the Abelian differentials. In this approach the role of the complicated 
hodograph system of equations is reduced to supply initial data for the system of first- 
order differential equations. With this simplified technique we are also able to derive 
several important general properties of critical processes Is — )■ Is±i such as the 
matching of their corresponding solutions of the Whitham equations or the (third) order 
of the phase transitions that they represent [3 [H [9l [TOj [23] . 

To illustrate our method we consider the Bleher-Eynard potential [7], [H] 
1 4 

V{x) = -x^ - -cx^ + {2c^ -l)x'^ + 8cx, -l<c<l. (10) 

This model was proposed as a concrete example of a third order phase transition I2 — ?■ Ii 
(merging of two cuts) at Tc = 1 + 4c^. However our numerical and theoretical analyses 
of the corresponding system of first order ordinary differential equations shows that for 
c 7^ the model also exhibits another third order phase transition Ji — t- I2 (birth of 
a cut) at a certain temperature Tc with < Tc < Tc. We also derive the asymptotic 
behaviour of the functions (3j{T) and d^F{T){n = 0, . . . , 3) near both critical points Tc 
and Tc. For the merging of two cuts we obtain the Bleher-Eynard asymptotic expressions 
for /3j{T) and calculate the magnitude of the jump of 9|.T(T) at Tc. For the birth of a 
cut we find the same logarithmic behaviours as in the Eynard models [9]. 

The layout of the paper is as follows. In section 2 we set up our notation and present 
the elements of the theory of Whitham hierarchies on hyperelliptic Riemann surfaces. 
Section 3 is devoted to the solutions of Whitham hierarchies underlying matrix models. 
In section 4 we discuss several general properties of Is — )■ Is±i phase transitions, whereas 
section 5 addresses the particular cases Ji ^ l2- We defer to Appendix A a self-contained 
proof of the relation between the free energy of the matrix model and the r-function 
of the hierarchy, to Appendix B some explicit expressions and asymptotic estimates of 
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elliptic integrals, and to Appendix C the derivation of asymptotic approximations to 
the solutions of the differential equations in neighborhoods of the critical points. 



2. Multi-cut matrix models and Whitham hierarchies 

We denote by w{z) the double- valued function 



w{z) 



2s 



^n^^-A), (11) 
\ i=i 

by Wi{z) the branch of w{z) with asymptotic behavior wi{z) ^ z^ as z ^ oo, and by 
the boundary value of Wi{z) on the real line from above. It has been shown |1] 
that the eigenvalue density p{x) takes the form 

Pix) = WiA^), xeJ, (12) 

Zm I 



h{z)={^] , (13) 



where 

h(7^ = ( 

and where {f{z))^ denotes the polynomial part of f{z) at z = oo. As a consequence 
of (jlj) the endpoints (3j of J must satisfy the equations |lj 

h{x)wi^+{x)dx = 0, j = 1, . . . , s — 1, (14) 
■ V'(z) 

.z^— ^dz = 0, j = 0,...,s-l, (15) 

where 7 is a large counter clockwise oriented loop which encircles J. Moreover, the 
normalization condition (E]) gives 

h{z)wi{z)dz = -AiriT. (16) 



7 



Equations (fT^ - (fT6|) give 2s conditions to determine the 2s unknowns f3i, . . . , (32s- 
However, it may happen that for fixed T and V{x) there exist different values of s 
for which the system f[T^ - (ITBi) has a solution /3i < ■ ■ ■ < /32s- If this is the case, the 
additional conditions p{x) > for all x E J, and ([5]) must be used to characterize 
the unique admissible solution of the problem. In terms of the functions h and Wi 
equation ([5]) can be written as [1] 

h{x')wi{x')dx' < 0, for x < (3i, (17) 



/ h{x')wi{x')dx' >0, for (32j < X < (32j+i, j = l,...s-l (18) 
h{x')wi{x')dx' > 0, for x > P2s- (19) 



132s 
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In addition there is an important a priori upper bound for the number s of cuts from 
a given potential V at any temperature T > 0. Equation f|T3|) shows that h{x) is a 
polynomial of degree 

degh = degV - s -1, (20) 

while from f|T4l) it is clear that h{x) has at least one zero on each interval /32j < x < P2j+i- 
Therefore deg h> s — 1 and we have the following bound: 

.<^. (21) 

Finally, the eigenvalue density p{x) is called regular if h{x) ^ for all x 6 J and the 
inequalities (fT7|) - (fT9|) are strict. Otherwise p{x) is called singular. 

A basic property of multi-cut matrix models in the large N limit is that (IHj)- 
(fT6|) . which determine the endpoints of the support of the eigenvalues, admit a natural 
interpretation in terms of Whitham hierarchies of quasiclassical integrable systems. 
In this section we describe the elements of the theory of Whitham hierarchies on 
hyperelliptic Riemann surfaces which are required for this interpretation. 

2.1. Hyperelliptic Riemann surfaces 

Let us denote by F the hyperelliptic Riemann surface associated with the curve 

2s 

w' = \{{z-P,). (22) 

i=l 

The two branches wi{z) and W2{z) = —wi{z) of the function w{z) characterize F as a 
double-sheeted covering of the extended complex plane: 

F = Fi U F2, T, = {Q = {wi{z),z) e F}. (23) 

The following coordinate systems constitute an atlas for F: 

kooi = for Q near oOj, i = 1,2 (24) 

kj = a/^--^, for Q near (0,/3j), (25) 
k = z, for Q near any Qo ^ {ooi, 002, (0,/3j)}. (26) 

where j = 1, . . . , 2 s. The homology basis {oj, of cycles in F is defined as shown 

in figure 1. 

The Whitham hierarchies describe certain deformations of the Riemann surface F 
induced by functions f3j = f3j(to,ti, . . .) in such a way that a suitable set {dflk}'kLo of 
Abelian differentials on F can be derived from a "potential" dS, i.e., dQk = dt^dS. 

We denote by dflo the third kind normalized Abelian differential whose only poles 
are at ooi and 002, and such that 

- + 0{z-^))dz, Q^ooi, 
dno{Q) = { }\ \ (27) 

- + o{z-^)jdz, g^CX)2, 
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where here and hereafter it is assumed tacitly that z = z{Q). For > 1 we denote by 
dQk the second kind normahzed Abehan differential whose only pole is at ooi, and such 
that 

dQkiQ) = {kz^-^ + 0{z-^))dz, Q ooi. (28) 
Note that the usual normalization conditions 

dfifc = 0, i = l,...,s-l, (29) 



determine the differentials df2fc uniquely [211 [22]. If is easy to see that 

dfifc = f I z^-i + dz, A; > 0, (30) 



.2 w[z^ _ 

where P^i^z^ are polynomials of the form 

40+2) (^'"'«^i(^))e + I]cfci^^ (31) 

' i=0 

and the coefficients C]~i are determined by the normalization conditions ( l29l) . We give 
explicitly the first polynomials for s = 1 and s = 2 which will be used in section 5. 

Case s = 1 For s = 1 we have 

P^z) = (40 + ^) {z'-' V(^- A) (2-/32))^, (32) 

so that the first few polynomials are 

Poiz) = 1, (33) 
Pi{z) = ^z-^iP, + (3,), (34) 

P2iz) = z'-^iP,+P2)z-^ (/3i - ^2?. (35) 
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Figure 2. Fundamental domain. 



Case s = 2 For s = 2 we have 

. , . , ^ k 

PkZ) 



\{z) = (^4o + 2] (^'~V(^ - - mz - mz - /34))^ + c,o, (36) 
where the coefficients Cfco are determined by the normahzation condition 

'''' ^'^"'^ dx = 0. (37) 



P2 



Using the elhptic integrals in Appendix B we find 

Poiz) = z-P, + iP,- Ps)^^^, (38) 

K{s) 

Pi{z) = -\{f3i+l32 + Pz + PA)z + \ (/3i /34 + /32 /Ss) 

+ ^(/34-/32)(/33-/3i)|||, (39) 

where K{s), E{s) and n(r, s) are the complete elliptic integrals of first, second and third 
class, respectively, and 

^ /33-/32 ^ (/34-/3i)(/33-/32) . . 

/34-/32' (/93-/9i)(/34-/32)' ^ ^ 

Our next step is to define Abelian integrals Qk associated with the preceding 
Abelian differentials in such a way that {dz ^k)dz = dQk- We first deal with the 
generic case {dQk}'kLi and then consider separately the special case of d^o- If we 
let all the basic cycles {ai,bi}^~l pass through some fixed point and cut T along these 
cycles we determine a simply-connected domain Fq, the so called fundamental domain 
of F (figure 2). Thus, by restricting our integration paths to Fq — {ooi} we can define 
single- valued Abelian integrals 

rQ 

^]fc(g)= / dn^ + Lk, k>i, (4i) 
J 002 
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where the 2;-independent term Lk is 



Lk = hm - / dfifc). (42) 

Q->ooi 



' 002 

Note that from (EHl) it follows that 



/ dnk = + 0{l), g-^ooi, k>l (43) 

J (X)2 

and therefore 

1^ Lfc + (^(Z ^), Q CX)2, 

In order to define a single-valued Abelian integral associated with df2o we fix a 
path Co in Fq connecting ooi with 002, and a complex number such that both regular 
points Qi = {wi{zo), zq) and Q2 = {w2{zo), zq) are in cq (figure 3). In this way we define 
the single- valued integral 

^o(Q) = r (^^0 + +Lo- log {z - zo), (45) 



where the integration path is in Tq — Cq, the constant term Lq is defined by 

Lo = lim (2 \ogz- r fdfio + ^^\] , (46) 

and log stands for the principal branch of the logarithm. Note that with this construction 

dz 

dfio + = Oiz-^) dz, Q -> 002. (47) 

Z- Zo 



Therefore the integral in fl45l) is convergent and path-independent in Tq — cq because 
all the poles of the integrand are in cq and consequently the sum of the corresponding 
residues vanishes. Furthermore, 



^rfdn„ + ^)=-^. »L„ = 0, (48) 

OZoJoooK Z-ZoJ Z-Zq dzo 
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and consequently flo{Q) is independent of zq. It follows easily that 

^'^^^ - \ -log Z + Lo + O{z-'), Q^oo,. ^^^^ 
2.2. Whitham hierarchies 

As we anticipated in the preceding section we consider deformations of the Riemann 
surface V induced by functions /3j = (3j{tQ,ti, . . .) such that 

dt^dnk = dt,dnj, j,k>o. (50) 

These equations are the Whitham equations for the differentials dflk and are the 
integrability conditions for the existence of a meromorphic differential dS satisfying 

dt^dS = d^k, k>0. (51) 

Hereafter we will use the vector notation j3 = (/5i, . . . ,(32s), and t = (^0,^1,^2, • • so 
that f3 = /3(t). 

Integrating ( ISTl) with respect to all the deformation parameters and enforcing the 
normalization conditions 

dS = 0, i = l,...,s-l, (52) 

we infer the existence of a unique normalized differential dS whose only poles are at cxoi 
and 002 with expansions of the form 



dS{t,Q) 



J2ktkZ^-^ + ^ + 0{z~^)]dz, Q^ooi, 

(53) 



z 

,k>l 



-^-^ + Oiz-^) ] dz, Q^ 002. 



Our next goal is to pass from this differential to its corresponding Abelian integral. 
Thus, from gl]), (|l2]), (iSl), (iE]) and ([50]) we obtain 

dt^Lk = dt,L„ j,k>0, (54) 

and therefore 

dt^ ilk = dt^ VLj, j, k>0, (55) 

which are the Whitham equations for the Abelian integrals f2fc(t,2;). Consequently, if 
we define by 

^(t, Q) = r (dS + to^^) + m - to hgiz - ^o) (56) 

J 002 \ Z — ZqJ 

the Abelian integral on To — Cq corresponding to the differential dS that satisfies the 
conditions (152|) - (!53|) . and if L is such that 

dt, L = Lk, k>0, (57) 
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then ( 15 ip implies at once that 

\Sit,z) = nk{t,z), k>0. (58) 

Moreover, by adjusting a constant term in L we deduce that the function S has the 
asymptotic expansions 

2_^tkZ +tQ\ogz - 2_^—-, Q -> ool, 

k>l i>l 

Vi -t^logZ + Q ^ 002. 



S{t,Q) 



I z' 

i>l 



The most basic object of the Whitham hierarchy is the quasiclassical r-function 
r(t) = expF(t). It is characterized [TlJ by the equations 

dF 

V^+l = z > 0, (60) 

where the functions Vi are the coefficients of the expansion of S given in f l59|) . We prove 
in Appendix A that 

^(t) = ^ (fv{z)dS{z) + Uov^ = ^J2t,v,+,{t), (61) 

where 7 is a large counter clockwise oriented loop which encircles J in Fi. Note that 
the calculation of the r-function requires knowledge of the vi coefficient, which in turn 
requires the construction of Abelian integrals and not just of Abelian differentials. In the 
next section we will see that F is precisely the free energy in the leading nonoscillatory 
factor in ([8]), and that Vi is the Lagrange multiplier in (j4]). 

3. Solutions of Whitham hierarchies underlying multi-cut matrix models 

The purpose of this section is to prove that given a family of matrix models with a fixed 
number s of cuts, then the differential 

dS=\ (v'{z) + (^^) ^w{z)^ dz = \ {V'{z) + h{z)w{z)) dz (62) 

is the normalized solution with asymptotic behaviour (!53|) of Whitham's equations for 
a finite number of deformation parameters to, ^1, • • • , i2p- Here to = —T and ti, . . . ,t2p 
are the coupling constants in the potential ([2]). 
First note that f|T^ simply means 



d^ = 0, i = l,...,s-l. (63) 



If we write dS in the form 



dS=l (v'{z) + ^M^) - f ^ ) w{z) ] dz, (64) 
2 V wi{z) \wi{z) ' 
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where (/(^))e = f{z) — {f{z))^, and recall that ^1(2;) ~ 2;'' as 2; — 00, then it is obvious 
that ([15]) and ([HD imply 

(^) w,{z) = -'^-^ + 0{z-'), z^oc, (65) 
\wi{z)Jq z 

or, equivalently, that dS has the asymptotic expansions 



dS{Q) 



+ ^ + j dz, asQ^ooi, 



(66) 

(^-^ + 0{z-^)^ dz, as Q ^ 002. 

Moreover, in terms of the local coordinate fcj near a point (0, Pi) we have that 

z = k^ + l3^^ dz = 2 kidki, (67) 

so that 

dS = 0{ki)dki, ki 0. (68) 
Hence the differential dS has poles at ooi and 002 only, and from ( 166|1 it is clear that 

2p 

dS = ^tfedfife. (69) 

A:=0 

Differentiating (IMj) with respect to (A; = 1, . . . , 2p), it follows that 

f {kz''-^ + 0{z-^)) dz, asQ->cx)i, 
dt.dSiQ) = I (70) 
I 0{z^'^)dz, as Q 002- 

Furthermore, f lBS]) shows that near (0, /3j) 

9t,d5 = 0{l)dkj, kj 0. (71) 
Hence, since the a^-periods of dt^dS vanish, dS satisfies the Whitham equations 

dt,dS = dnk, fc = l,...,2p. (72) 

The case k = can be proved similarly. 

As a consequence of (jHSD and fpU]) we obtain an important identity which for later 
reference we state making explicit the dependence of P, h and w on the endpoints f3: 

J2 tkPkiz, (3) = h{z, f3)w{z, Pf. (73) 



fc>0 
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3.1. Hodograph equations 

If we set z = f3i in ( ffSl) . we find that the solutions of the Whitham equations ( ISTI) 
associated with multi-cut matrix models satisfy the hodograph type equations 

2p 

Y,tkPkW^,f3)=0, t = l,...,2s. (74) 

fc=0 

These equations can be stated in a geometric formulation fl2[ [13] if we recall that in 



terms of the local coordinate ki{t) = \/ z — (3i{t) near (0, /3i), we have that for all /c > 
dfi, = f ^^^^^"^^ + Oih) I dh, h ^ 0. (75) 

Hence, equations ( l69l) and (1741) imply that 

d5|(o,ft) = 0, t = l,...,2s. (76) 

Thus the solutions of fl74|) determine zeros of the differential dS. Moreover, equating 
the principal parts of dt^dQk and of dt^dfli at the poles (0, /Si) and taking into account 
that dt„dki = {dt„/3i/2kf)dki, we find that (17i|) represent the hodograph transform [26] 
for a compatible set of diagonal hydrodynamical systems given by 

These are the Whitham equations in hydrodynamic form. 



dtA = \:)7 JA A. k > 0. (77) 



3.1.1. Example The simplest non-trivial example that we may consider is the quartic 
even potential [6]: 

V{z) = \z' - z\ (78) 

Its associated hodograph equations ( !7i|l are 

^ P4(A, /3) - P2(A, /3) - T Po(A, /3) = 0, « = 1, . . . , 2s. (79) 

According to the bound ( l2Ti) . as T changes there are only two possibilities for the number 
of cuts: s = 1 or s = 2. For s = 1 the hodograph equations (!7il) read 

-256T + 35/3f-5/3^-20 /^j^ /32 -6 /^^^ /^s^ -4 /3i -96 /3i2 + 32 /32^ + 64 /3i /^a = 0, (80) 

-256T + 35/32^-5/3f-20 Pi (3i-Q /3i /3f -4 /32 ^! -9Q /3i + 32 ^l + 6A /32 /3i = 0, (81) 

and the function h{z) is given by 

h{z) = z' + ^{f3, + (32) + ^{3f3l + 3 Pi + 2 /3i /^a) - 2. (82) 

Imposing the symmetry condition = — the hodograph equations reduce to the 
single equation 

2T + - = 0, (83) 
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which leads to the solution 



1 + Vl + 3T, (32 



1 + Vl + 3T, 



4) 



' '^^ 

which is real and satisfies Pi < P2 for all T > 0. However, the corresponding function 
h{z) takes the form 



h(z) 



+ - (1 + V1 + 3T) - 2, 



^5) 



and, obviously, is strictly positive on the real axis only for T > 1. Furthermore, we have 

' -\x^-/3f\'/\ forx</3i. 



= <^ i - for (3i<x<l32 



(86) 



^ Ix^ -/^i^/^ for (32 < X. 

Hence for T > 1 the functions (|8^ are the endpoints of the support of eigenvalues, since 
p{x) > for all X G {f3i, (32) and the inequalities (fT7|) -(fT9l) are strictly satisfied. 

If for s = 2 we impose the symmetry relations = — (3^ = —(32, the terms Cko 
with even k in (1361) vanish and the hodograph equations (1791) reduce to the system 



/3i ( 2T - l(3t - (3i + 1(31 + (3! + \m 
/32 ( 2T + \pl + Pi - ^-p', - Pi + 1/3,2/32 



0, 
0, 



which supplies the following solution of f[74j) : 
/3i = -a/2(1 + Vf), P2 = - 



2(1 -Vt), 



Ps=^2{l-Vf), p, = ^2{i + Vt). 

Moreover, in this case h{z) = z and 

f I (x^ - /32) (x^ - /32) 1 1/2^ for X < /3i or X > /34 



(87) 
(88) 

(89) 
(90) 



wi,+ (x) 



-I |(x2 - Pl){x^ - /32)r/2^ for /3i < X < /32, 
-|(x2 - Pl){x' - Pl)\"\ for /32 < X < /33, 



(91) 



[ z |(x2 - Pl){x^ - /32)|^/2^ for P,<x< P,. 

Hence for < T < 1 the functions flS^ - flUU]) are the endpoints of the support of 
eigenvalues, since p(x) > for all x G {Pi, P2) U {Pzi Pa) and the inequalities (fT71) - (fT9ll 
are strictly satisfied. Therefore there is a merging of two cuts at T = 1. Note finally 
that the functions {Pi, P2) of the s = 1 sector have a C^-smooth matching with the 
functions {Pi, P4) of the s = 2 sector at T = 1. Indeed, for s = 1 we have 



Pi = -2--{T-l) + 1(T - If + 0{{T - If). 
4 io 



(92) 
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Figure 4. Evolution of tlic cndpoints for tlie quartic even potential. 



while for s = 2 

/3i = -2 

Notice that the matching is but not C^. 



\{T-l) + ^^{T-lf + 0{{T-lf). 



(93) 



3.2. An alternative approach to determine the eigenvalue support 

Except for the simplest examples, the system of hodograph equations (!74|) is extremely 
difficult to solve, even by numerical methods. Our strategy to deal with this difficulty 
is to take advantage of a particular feature of the solutions of the Whitham hierarchies 
associated with matrix models: the functions /3j(t) satisfy a system of first-order ordinary 
differential equations which will provide not only a practical numerical method but also 
analytical estimates for studying the asymptotic behaviour near critical temperatures. 
From ([62]) and (El) it follows that 



where n 

we get 



0, 



5.„log(M.,/3K(.,/3))=2.^^^^^^^^^^^^^„ 
. . . , 2p. Then, equating the residues of both sides at z 



4P„(A,/3) 



(94) 
,2s) 

(95) 



MA,/3) n,y.(A.-/3.)' 

By construction this is a compatible system of first-order ordinary differential equations 
which is satisfied by the endpoints functions. Therefore, if we can find a solution /3(to) 
of the hodograph equations (1741) at a certain value to of the deformation parameters, 
then this /3(to) will supply a set of initial data for solving the system ( l95l) and determine 
completely /3(t) in a region of constant s. 

We stress that the first order ordinary differential equations ( 195|) in fact imply the 
full system of hydrodynamic Whitham partial differential equations (177|) . Furthermore, 
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if we denote by Xi, . . . , Xs-i the roots of Po{z, (3), then the n = member of reduces 
to 

w,/3)n,v.(A-/3,)' ^ ^ 

which are the Bleher-Eynard relations [A. 13] in [7]. 
3.3. Relation between the r -function and the free energy 

We now estabhsh the relation between the r-function and the nonoscillatory factor of 
the partition function of the matrix model in the large N limit. From (1691) and ( 1721) we 
deduce that 

J2tkdudnk = 0. (97) 

As a consequence the functions Lk defined in fH2|) and fH6|) satisfy 

dt„L{t) = L^, n = 0,...,2p, (98) 

where 

2p 

L(t) = ^tfcLfc(t). (99) 

fc=0 

Thus the Abelian integral fl56|) corresponding to the differential fl62l) is 

2p 

5(t,Q) = ^4fifc(t,Q), (100) 

fc=0 

and has the asymptotic behaviour given by fl5^ with 

t;i(t) = L(t). (101) 
In particular we get that Vi satisfies 

dt„ Vi = L^, n = 0, . . . , 2p. (102) 

To relate the function 5* to the eigenvalue density p{x) of the matrix model we prove 
that the restrictions Si of the function S to the sheets {i = 1, 2) of the Riemann surface 
can be written in the form 

5i(t, z) = V{t, z)+to j p{x) log(z - x)dx, (103) 

S2it, z) = vi{t) - to j^p{x)log{z - x)dx. (104) 
We begin by noticing that the derivatives of Si with respect to z are given by 

s;w = i(v'w + (^)^ (105) 

Si{z) = V'(z)-SHz). (106) 
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Thus, the boundary values of S[{z) on the real axis satisfy 



Hence 



Wi{x), for X G J 

S[ix + iO) - S[ix - iO) = { ^"^1^^^^® (107) 

0, for X G (— oo, oo) — J. 

S[{z) = V'{z) + to [ -^dx, (108) 
Jj z-x 

S'^{z) = -to j^-^^dx. (109) 

Therefore, taking into account the normalization of p{x) and fl59|) . our equations (11031) 
and (11041) follow. Moreover since the function 5* is continuous in Fq — Cq (cf. figure 3) 
we have 

Si{x + iO) = S2{x-i0), for all xG J. (110) 
Thus, using (fT03|) . (fTOlj) and 

log(x — ?/ + iO) + log(x — y — iO) = 2 log |x — yl, (xj^y), (HI) 

we obtain 

v^{t) = Vix) + 2to j^p{y)\og\x-y\dy, x G J. (112) 

Therefore 



t>i(t) = / t>i(t)p(x)dx = / y(x)p(x)dx + 2to / / p(a;)p(y) log |x — ?/|dxdy. (113) 
J J J J J J J J 



Moreover, (1121) and ( ]62l) imply that 

1 
2 



Viz)dSiz) = \(b V{z)h{z)w{z)dz (114) 



V {x)h{x)w ij^{x)dx (115) 

= 2Tiito J V{x)p{x)dx, (116) 
and we are led to the main result of this section, namely 
F(t) = log r(t) = ^ V{z)dS{z) + 

= — T / V^(x)p(x)dx + / / log |x — ?/|p(x)p(y)dxd?/, (117) 
J J J J J J 

where in the last equation we have substituted to by — T. 
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4. Phase transitions in multi-cut matrix models 



For a given fixed potential let us denote by Is the set of values T > such that the 
eigenvalue support J has exactly s cuts and the corresponding eigenvalue density is 
regular. These sets correspond to the different phases of the matrix model. Common 
boundary points of two sets Ig and Is-i mark phase transitions in which the number 
of cuts changes in one unit. Our next goal is to match pairs of solutions of Whitham 
hierarchies to characterize critical processes. 

Hereafter we will use a superscript '(s)' to distinguish objects corresponding to 
different multi-cut cases. Using our new, more explicit notation, we recall that for each 
T E Is the endpoints of the eigenvalue support 

Pi'\T)<(3i'\T)<...<(3i:\T), (118) 

satisfy the hodograph equations 

2p 

rp(^)(/3f\/3(^)) -^t,Pf (/3f),/3W) = 0, . = 1,...,2.. (119) 

k=l 

As we mentioned in the preceding section it is extremely difficult to characterize 
(3'''^\T) from this system. Thus, we will take advantage of the n = member of the 
system (195|) 

A'^ = < \ ' M rv^ fc = l,...,2., (120) 

/.W(/3(^),/3(^))n.^,(/3i^)-/3fV 

where the dot denotes derivative with respect to T. 



Phase transitions with Ig — )■ Is±i 

In this section we will use matched solutions of Whitham hierarchies to describe a type 
of critical processes in which two consecutive endpoints of the eigenvalue 

support coalesce at a critical temperature T = Tc, and such that: 

(i) The remaining components of /J^"-* and of ^ match continuously with the 
corresponding limits of f3^'^~^^ and /3*" 

(ii) The first and second derivatives OtF and OttF of the free energy are continuous 
at T = Tf.. 

If we denote by and (3^ the limits of the endpoints and their T-derivatives as 
T ^ Tc, these conditions mean that 

^ = p%., {AT\ptu'^) = \ (121) 

I (/3S+2,/3S+2), l<k<2s-2. 
As an example of this notation, the phase transition in figure 4 corresponds to the case 
s = 2, / = 2 in (11211) . where the red labels correspond to the '(2)' phase, the blue labels 
to the '(1)' phase, (i'^'' merges with [if \ pf'' matches and (if^ matches 
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As we will see, our method to characterize these processes is based on the important 
identity (I73|) and the following property of the polynomials Pk (see Theorem 3.5 of [T9]): 
given (s > 2) such that for some / = 1, . . . , 2s — 1 

A = A+i = /3, (122) 
then the following identities hold 

Pt\z,(3^'y) = iz-P)Pt'\z,f3^'-'^), k>0, (123) 

where 




(124) 

I <k<2s-2. 

We give some simple examples of identities f ll23p which can be directly verified and that 
will be used later. 

4.I.I. Examples Consider the polynomials Plf\z, (3^^^), {k = 0,1) with = 
{(3i, (3, (3, (34). From the expressions (l38i) and (!39l) . taking into account that (32 = 
implies r = s = 0, and that 



E(0) = ir(0) = n(0,0) = ^, (125) 
we deduce at once the identities 

Pl,'\z,(3^'^) = z-^ = {z-P)Pi'\z,(3^'^), (126) 

pl'\z,f3^'^) = iz-(3)(^^z-^i(3, + (3,)^={z-(3)Pi'\z,f3^'^). (127) 

Likewise, consider the polynomial Pq'^\z, (3^'^'') with /S*-^^ = (A, A, /?)■ Now 
= Pi implies r = s = 1, and from the asymptotic approximation fl27ip of Appendix C 
it follows that 

Hence, (!38l) implies 

p(^)(^,/3(2)) = ^_/3 = (^_/3)P«(;,,/3(i)). (129) 

4- 1-2. -matching of eigenvalue supports For the critical processes that we are 
considering the limit values (3^^\ {q = s, s — 1) of the endpoints of the eigenvalue support 
satisfy the hodograph systems 

2p 

TA'\(i'^$.(i^^) -Y.^k Pt\PS^f3i'^) = 0, ^ = 1, . . . , 2g. (130) 

k=l 

Since /3^^^ verifies the hypothesis f ll22p then from (11230 we deduce that 

Pt\^,f3i'^) = {z-(3)Pt'\z,(3t''>). (131) 
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As a consequence we have 

Pt\Pi% f^i'^) = - /5) Pt'^iPS^f^t'^)^ ^ = 1, . . . , 2.. (132) 

Hence, the hodograph equations (11301) for g = s and i = 1,1 + 1 are trivially satisfied, 
while the remaining equations reduce to the system (I130p for q = s — 1. Thus, the 
hodograph systems are compatible with a continuous matching of the endpoints of the 
eigenvalue support, 

{k, l<k<l, 
(133) 
k-2, l + 2<k<2s. 

We notice also that fl73|) . and (I13ip imply 

/»'-'(../3[-')^ """;':f"' . (134) 

and therefore the continuity of the eigenvalue densities at the critical point follows 

pW(x,/3(^)) = p(-^)(x,/3^^)). (135) 
Furthermore, from (I120p . (11320 and (I134p we deduce that 



Ms) _ ^^0 \Pc,k^Pc 

Pc,k 



= /3(ji), A;^/,/ + l, (136) 
which shows the property for the matching of the eigenvalue supports. 

4-1.3. Order of the phase transitions We will now prove that these critical processes 
are, generically, third-order phase transitions. 

Consider the behaviour of the derivatives of the free energy (I117P with respect to 
T as T — )■ Tc ±0. From the expression of the free energy (I117p and from (I135p it is clear 
that F is continuous at T^. Now fl60l) and (I112p imply 

dF{T) 



dT 



— f 1 = —V{x) + 2 T J p{y) log \x — y\ dy, x E J, 



(137) 



-drvi = Lo= lim (2 \ogz- ( dQo + ] ] . (138) 



Therefore OtF is also continuous at and from (130|) and (I123P it follows that 
92fW(T,) d^F^'"'^{T,) _ f^' [Pi'-^\z,l3'f-^^) Pif\z,f5^''^ ' 



dT^ dT^ J 00, \wi^~^){z,(3i'-^'^) w^^){z,(3['^) 

which shows that OttF is continuous at 
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Finally, as a consequence of 



s^F r'ad.„ r'j2'-^^, (140, 



and using again fl30|) and (I123p . 

^(-,/3i^^) = ^^(-,/3r^ jVM + 1. (141) 
Therefore, the jump of the third derivative is 



T^T. \^ Q^is) Q^is)^ 

As we will show with an explicit calculation in section 5, this jump is in general different 
from zero and, consequently, these processes represent third-order phase transitions. 



5. Phase transitions Ji ^ I2 



In this section we first particularize our previous analysis to the case in which Tc is a 
common boundary point of two regions Ji and I2 and then apply these results to the 
matrix model corresponding to the Bleher-Eynard potential [?]. 

For s = 1, fl32l) implies that the hodograph equations flllQp can be written in the 

form 



2T 



^-V\z) ] M;i'^(z,/3«^ 



0, ^ = 1,2, 



(1) 



where 



«;f)(z,/3«) = 7(z-/3«)(z-/3«). 
Moreover, fll20p reduces to 



(143) 

(144) 

(145) 
(146) 



Similarly, from ( 136|) it follows that the hodograph equations flllOp for s = 2 can be 
written 
'2T 



(1) 



+ f/(T,/3(2)) = 0, 2 = 1,2,3,4, 



where 



wf\z,(3^^) = ^J{z-P?){z-pf\z-pf\z-pf^), 



(147) 



(148) 
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and 



Although these hodograph equations are extremely involved, the corresponding ordinary 
differential equations fll20p for /3*^^^ are much simpler. Indeed: 



where 



/o(2) _ n(2) 



and 



'"(/jf-/jP)(/jf-/jf)' ' ' 

C(/3l=') = - (0f ' - Pi>) SM. (153) 

As we have seen in the preceding subsection the second derivative OttF of the free 
energy (the specific heat) 

lim (2 log z - r (dQo + ] , (154) 



is continuous at T = Tc. Then in a transition Ji ^ /2 it can be calculated from the 
limit value of OttF^^^- In this way, integrating from 002 to ooi along to the path shown 
in figure 5 and taking into account that 

dfii^^ = , =, (155) 

^J{Z-P[''^){z-Pi''>) 

we get 

^^2,0.^. (15.) 

5.1. Merging and birth of cuts in the Bleher-Eynard model 
The Bleher-Eynard potential [7] 

V{x) ='^x^ -^cx^ + {2c^ -l)x^ + 9,cx, -l<c<l, (157) 



Phase transitions in multi-cut matrix models and Whitham hierarchies 



22 




Figure 5. Integration patli from 002 to ooi. 



defines a matrix model exhibiting interesting critical processes. According to the upper 
bound ( 12T|) there are two possible phases s = 1,2 only. In these cases it follows easily 
that 

2 



= \z-2c + 



+ Pi'' 



-2/3i'^/3^'^-c(/3r+/3r)-2, 



and 



(158) 



(159) 



Merging of two cuts We first look for a particular solution of the hodograph equations. 
A direct calculation shows that 

( (t " V^'(2)^ w']'\z, /3(^))^ =2T-z^ + 4cz{z^ - 4) - 4(0^ - l)z^ + Sc" - 2, (160) 

and as a consequence it is easy to see that at Tc = l + 4c^ the hodograph equations (11431) 
are satisfied by /B^^^ = (/3c\\/3c2^) = (~2,2). Moreover, we have 

/i(i)(^,/3(i)) = (^-2c)2, (161) 

so that the eigenvalue density vanishes at x = 2c G J. Hence Tc is a a common boundary 
point of Ji and I2, and the endpoints (/^g , (3^ ) of the eigenvalue support for s = 2 tend 
to /3 = 2 c as T — 7- Tc in I2. Thus 



(2) _ 0(2) 



J (2) 
c,4 



and 



r(/3f) = .(/3f)=0, /.(2)(z,/3 



(2)^ 



with -2 < /3 = 2 c < 2. 
From fll45p we obtain 



-1 



2, 



z-(3. 



1 



(162) 

(163) 



(/3 + 2)2' ^-^'^ ^^'2 (/3-2)2- 



(164) 
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If we write 

/i«(z,/3«)= ^^-2c + ^i^^^^j +a{/3^'^), (165) 

where 

= ^ + - 1 - c + - 2, (166) 

and use the expressions (fT64l) for the velocities /^^^^^ (« = 1, 2) we get 



^"(/3'"(r)) 



This imphes that for T shghtly to the right of Tc the function /;,*^^)(x, /3*'^^(T)) is strictly 
positive for all real x and T is in Ji. Consequently, if T slightly to the left of then T 
is in I2. Therefore, the critical process at Tc represents a merging of two cuts. 

From the expressions ffTMD for the velocities /3^7 it follows that as T - 

/3f)(T.-t)~-2 + ^^^, /3f(^c-t)~2-^^4^, (168) 

whereas the expressions (nm for I3r and become undetermined because — 
C{j3^'^^), h{(if'\(3) and /32^^ — /33^'' vanish. However, we show in Appendix C 
(equations fl247l) -fl257D) that 

C(/3(^))~/3f + ^(/3f -/3f), (169) 
and therefore fll5Up leads to 

~ft(/3f'./3'^')(/jf-/3p')(/jf-/jf)' 

which imply the following asymptotic behaviours [7j: 

pf\Tc-t)^P-^^!l= pi'\n-t)^/3 + ^^ (172) 

Analysis of the third- order phase transition at T = Tc We have that 
d^F p^Sdfio d fPo{z,/3) 



dT^ J 00, dT J^^ dT \w{z,(3) 

and it is clear that 



d^, (173) 



(174) 
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If we now use the Taylor expansion at /3^^'' = {f3^\ pf\ f^^'^) 
7(z-/3f)(z-/3f) = (^-/32). 



(2) 



.-/3,-^^^-i^^^||^ + 0((;3f -/3f)3), (175) 
^ 8(2; - Pi 



^2 



then, taking into account fl244p . we have that as T —)■ Tc 

^J {z - (5f^){z - pf^) -{z- C{(3^'^) ~ i (--1^ + (/33 - /32)^ (176) 

Hence it follows from f ll72p that 

Po^^^(^,/3«) Po^'^(^,/3(^)) 2(r-T,) / 1 P 



w;«(^,/3W) ti;(2)(z,/3(2)) (4-/32)^1^ V(^-/3)' (4 - /32)(^ - /3) 
Therefore integrating from 002 to ooi along the path shown in figure 5 we get 

4 f f°° dx /3 dx 



(177) 



4-/32 yj^ _ y^^34 4 - /32 72 (a; _ /3) 
4 



:i78) 



(4-/32)2- 

This equation shows that the third-order derivative of the free energy jumps at the 
critical point 



Numerical calculations In section 3.1.1 we have solved directly the simple hodograph 
equations (I80l)- (!8T1) for the symmetric quartic potential (|78il in the s = 1 region and, 
by imposing in advance the symmetry relations between the roots, we have solved also 
the hodograph equations (l87 I) -(l88 l) for the same potential in the s = 2 region. The 
symmetric quartic potential ( ITHj) is the particular case c = of the Bleher-Eynard 
potential f ll57p . but nonsymmetric (i.e., c 7^ 0) cases of this latter potential are not 
amenable to the same direct procedure. The corresponding hodograph equations in the 
s = 1 region are still purely algebraic, and with the help of a symbolic manipulation 
program the solutions can be found implicitly among the roots a two- variable polynomial 
equation of degree 12. The hodograph equations in the s = 2 region are already 
nonalgebraic (the roots appear in the arguments of elliptic functions) and even a 
numerical solution and identification of the roots seems to be extremely difficult. 
However, we have shown that the roots also satisfy the system of first order ordinary 
differential equations f ll20p . which is readily amenable to numerical integration. For 
concreteness we will outline the procedure in the particular case c = 1/2, which turns 
out to be typical, and omit the superscripts (cf. figure 6). 

We take as our starting point the critical temperature = 2 (marked by the 
rightmost dashed vertical line in figure 6), for which we have already discussed the 
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solution of the hodograph equations. We have also shown that sufficiently small values 
of T to the right of belong to an s = 1 interval, and that at T = Tc the system f ll20p 
for s = 1 is regular (the denominators do not vanish). Therefore for T > Tc we can 
integrate the system fll20p directly with the initial condition (3i(Tc) = —2, f32(Tc) = 2. In 
this interval the system is numerically well behaved and we obtain the two curves shown 
in figure 6. We have checked that within our numerical precision these solutions coincide 
with two arcs of the algebraic curve mentioned in the previous paragraph. Note also 
that the mere success of the numerical integration does not guarantee that the system 
remains in an s = 1 region and we have to study later the behavior of the density of 
states. 

Likewise, sufficiently small values of T to the left of Tc belong to an s = 2 interval, 
but at T = Tc the system ( 1120p for s = 2 is singular and we cannot start our integration 
to the left simply with the initial condition /3i(Tc) = —2, f32(Tc) = /33(Tc) = 1, 
l3i{Tc) = 2. We circumvent this problem by taking as our initial point a value Tq = T^ — t 
with t sufficiently small, and using the asymptotic behaviors f ll68p and f ll72p as initial 
values. As the integration proceeds to the left we find that the solutions /53(T) and P^iT) 
coalesce at a second critical temperature Tc ~ 1.845 097, where we find numerically a 
birth of a cut. This second critical temperature is marked by the leftmost vertical dashed 
line in figure 6. Incidentally, the computational effort required by the adaptive numerical 
integrator increases as the integration approaches the Tc, where the system of differential 
equations is again singular and the arguments of the ensuing elliptic integrals tend to 
one. We have found that we can achieve higher precision if at an intermediate point, 
say T = 1.85, we rewrite the system of differential equations using the relations [281 129] 

K{s) = ^^k(^\, |arg(l-s)| <7r, (179) 



n(r,s) = - ] n (^,^) , |arg(l-r)| <7r,|arg(l-s)| <7r, (180) 

(l-r)vl-s \r - 1 s-lj 

to map the branch points of the elliptic integrals at 1 to — oo. 

Finally, from Tc to T = we integrate again the nonsingular system of differential 
equations fll20p for s = 1 taking as initial conditions the values we have obtained in 
the previous step for /3i(Tc) and /32(Tc), and get the two curves shown in figure 6 that 
coalesce at T = 0. As a double check of this last step we point out that, within our 
numerical precision, the solutions again coincide with two arcs of the of the degree 
12 algebraic curve mentioned previously. Moreover, the point /3i(0) = /32(0) can be 
calculated exactly, because it is the position /3min of the minimum of the potential V"(x), 
i.e. the real root of the cubic equation 

- 4c/3^i, + 2(2c^ - l)/3^in + 8c = 0, (181) 

or, using the Cardano formula, 

4 2^/3 . N 1/3 4c2 + 6 , , 

/3min = 7TC+ — (-4c=^-36c + 3y6v^) + TIT^' (182) 

21/33 (^_4c3 _ 36c + 3^/6 VAj 
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Figure 6. Evolution of tlic cndpoints for tlie Bletier-Eynard potential with c = 1/2. 



where A = 4 c"^ + 22 



1. In particular, for our case c = 1/2 we find 



'mill 




- - -\/37-3VTu] 
\/37 - sVTu J 



1.269 53 



(183) 



in agreement with our numerical results. 

We next use ( fT2l) and (fT3l) to calculate the density p{x) in each of the three 
subintervals < T < Tc {s = 1), Tc < T < Tc {s = 2), and Tc < T < Tmax (■5 = 1), 
and verify that p{x) > for all x G J and that p{x) is normalized. As a sample of 
our numerical results, in figure 7 we plot the density as a function of the position x for 
the first critical temperature Tc ~ 1.845 097, for a temperature T = 1.9 in the s = 2 
interval, for the second critical temperature = 2 and for a higher temperature in the 
second s = 1 interval T = 3. Within our numerical precision all these densities are 
normalized. The transition from the first to the second plot illustrates clearly the birth 
of a cut; as the temperature increases both intervals lengthen until at = 2 they join 
in the third plot, illustrating the merging of the two cuts; finally, in the last plot we see 
the normalized, one interval density after the two cuts have merged. 

Finally, in figure 8 we illustrate the sufficient conditions ( fTTl) . ( fTSl) and ( fT9l) . The 
figure corresponds to the second plot of figure 7, i.e., to a temperature T = 1.9, for 
which the solutions are (3i ^ -1.989, /^s ^ 0.646, /Jg ^ 1.431 and /S^ ^ 1.870. Since this 
temperature belongs to the s = 2 interval we have the three conditions. 





(184) 



(185) 



for X > Pi 



(186) 



which correspond respectively to the three plots. 
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Figure 7. Density p{x) as a function of the position x for the Bleher-Eynard potential 
with c = 1/2 at four temperatures. The first two plots illustrate the birth of a cut at 
Tc w 1.845 097; the last two plots illustrate the merging of two cuts at Tc — 2. 




1.9 2.0 2.1 2.2 2.3 2.4 2.5 



Figure 8. Illustration of the sufficient conditions (|184p . (|185l) and (|186l) for the Bleher- 
Eynard potential with c = 1/2 in the s — 2 region (T = 1.9). The approximate values 
of the endpoints are /3i « -1.989, ^2 « 0.646, ^3 « 1.431 and /34 ~ 1.870, and the 
plots show the respective integrals as a function of x. 



Birth of a cut We have seen that temperatures shghtly below belong to a two- 
cut region 12- In the particular case c = where Bleher-Eynard potential reduces to 
the quartic even potential studied in 3.1.1, in which this two-cut region is the whole 
interval < T < Tc. However, for c = 1/2 the numerical integration of the differential 
equations f ll20p has revealed a birth of a cut at a certain temperature Tc between and 
Tc. According to the variational equations @ and ([5]), as T — )• the closure of the 
eigenvalue support reduces to the set of absolute minima of the potential V{x). But, 
except for c = 0, this set is given by the single point /3min, and then a birth of a cut 
must happen at a certain < Tc < Tc for all c 7^ 0. 

We denote the endpoints of the support at the critical temperature T = T^, by /3c,i, 
Pc,2, and = = P- Note that 

ripT) = sipT) = 1, (187) 
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and 

(z, ) = ^ + ^ 0S + fS) + 3^ - 4c. 
We prove in Appendix C that as T Tc 



n(r,s) 2^(/3-/3,,i)(/3-/3c,2) 



tanh ^ 



c,2 



As a consequence, 
Thus ffTiSj) leads to 



c,l 



In order to characterize the behaviour of (3^ and (3l we use fllSOp . which 
imphes 

^(2) _ Ajl- U{r,s)/Kis)) 



M2)_ AIl{r,s)/Kis) 



Using f lT89D it is easy to see that ffT93D and ffT94D imply that 
where 



7 



'4/3 + + /3,,2 - 8c) ^(/3-/3,,i)(/3-/3,,2) 
In this way, solving fll95p and f ll96p yields P 



/3f(T, + t)~/3-(-27 



logt 



1/2 



/3f(T. + t)~/3+(-27j^ 
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Analysis of the third-order phase transition at T = T^. In this case the difference 
between the third-order derivatives of the free energy above and below the critical 
temperature diverges [9]. Indeed, proceeding as in the analysis of the merging of two 
cuts we have 



OOl 



002 



ddno 

or 



OOl 



d fPo{z,(3) 



002 



dT \w{z,f3) 



dz. 



(200) 



p(^)(z,/3«) P;,''{z,(3^-^>) 



(2), 



7(z-/3f)(^-/3f) -{z- C{(3^^) 

(2) (2) (2) f2) 

We now use the following Taylor expansion at ' = (/3{ ,132 , (3^ , ] 



{z-/3P){z-/3r>) = iz-/3r) 



(2)^ 



j(2)^ 



z-a 



(2) 



5(2) n{2) 



Furthermore, (11891) implies that 



2V(/3-/3,0(^-^,,) 

0-/3 + , .^(2) ^(2), tanh 

log(/3l - /^s ) 

Hence, we deduce that 

P(^)(.,^W) p(^)(z,/3(^)) 7 



13 -(3, 



c,2 



\ /3 - /3c,i 



wW(z,/3«) w;(2)(^^^{2)) log(r-re)^ 



where 



4V(/3-/3c,i)(/3-/3e,2) 



7 = --^ ' ^ ' ^ = tanh 

V(^-/3c,i)(^-/3c,2) 

Therefore we find that the integrand of 
diverges as (T - TJ-^ (log(T - f,))-^. 



(201) 



(202) 



(203) 



(204) 



(205) 



(206) 
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Appendix A: The r-function of the Whitham hierarchy 

In this appendix we prove that the free energy F defined in (1611) satisfies (jSO])- The 
derivative of F with respect to tn for n > 1 is 

= ^ii'^"^^^'^ + ^J^V{z)dtMS{z)) + ^jdt^v,. (207) 

We analyze separately each of the three terms of this expression. The first one is easily 
evaluated using f l59|) : 

The calculation of the remaining terms is more involved. Using (15T|) we can write 
the second term in the form 

^ j( V{z)dt^ {dS{z)) = ^ V{z)d ^n = ^^Y.^tmj^ z^dn^. (209) 
Now, taking into account fl28l) and fj44|) . 



= - f (l^n)ed(f^m)e = - f ^nd{Q.r, 

J 'y J y 



= - (j)nndnm+ f^nd{nm)e- (210) 

The first integral in the last expression can be transformed into an integral along the 
boundary STq of the fundamental domain (see figure 9): 

VLr,dVLra = -i ^n<^^m- (211) 
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But using the normalization properties and the bihnear relations for Abelian differentials 
dfifc [211 125] we have that 

(f ^nda^ = ^( (f d^n ■ i dVlm - (f dfi^ ■ (f d^n) = 0. (212) 

JSFo i—i J ai Jbi Jai Jbi 

Moreover, it follows easily that 

(fn^d (nje = / ^" d {nje = / ^" d n^. (213) 

1/7 

Hence, (I210p implies that 

Z'^drin = I z'^dVtm, (214) 
J 'y 

and using (|69l) and fl209p we finally obtain for the second term 

^ ^ vi.)ru<i Si.)) = ^ d Si.) - ^ I n„ d n„ 

= ^-^/n„dn„. (215) 

To calculate the third and last term in (120 7p note that from (11131) . (jH]) and (!27|l it 
follows that 

dt^Vi = Ln = -77^-7 / ^ndVto, (216) 

27ri 

where 7' is a large counter clockwise oriented loop which encircles J in Now (I208p . 
f l2T5|l and f l2T6|) imply 

= - ^ / l^ndl^o, (217) 
47ri 

and integrating along 5Vq (see figure 10) and using again (!29|) . we have 

f2„dfio — ^ f ^nd^o 
7+7' "''5ro 



^ r i" dfi„ ■ / d^]o - / dfio ■ / dfi™") = 0. (218) 

j—]^ \ J a; Jbi Jai Jbi J 



Therefore we conclude that 

dt^F = Vn+i, n>l. (219) 
Consider now the derivative of F with respect to to. From f l5T]) and f lll3p we have 

dt,F = ^ V{z)dno + J + |Lo. (220) 

Moreover, following the same procedure as in the proof of (12191) . 

n>l A n>l A' ^ „>i 
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Figure 10. Paths 7 and 7'. 



and therefore, using fllOip we obtain 

1 Vi to 



ra=l 



which completes the proof of (160|) . 
Appendix B: Elliptic integrals 

The complete elliptic integrals of first, second and third kind are given by 

dt 



K{s) = 
E{s) = 
n(r, s) 



vT^Vi-st^' 
^ VI - st^ 



■dt, 



dt 



Vr^(l-rt2)Vl-st2' 



respectively. In terms of these integrals, 

dx 



2K{s) 



J 02 



^{x 




xdx 


-mx 


-Pa) 






-/32)ix 
x'^dx 


-/33){x 


-Pa) 



2{^^K{s) + {P^-P,)Tl{r,s)) 
V(/33-/3i)(/54-/52) ' 



/-5 «TV (2/34 - Wa - m^A - P2))K{s) 

V IPs — Pi)\Pa — P2) 

- rni^A - /32)E{s) - (A + /32 + /33 + /34)(/34 - /3s)U{r, 
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In our study of the merging of two cuts we need the Taylor series of n(r, s)/K{s) 
near (r, s) = (0,0). From the integral representations (12231) and (12251) it follows 
immediately that |27] 

TT / s 9s^ \ 
Kis) = -1 + - + — + ... , s^O, (229) 



2 V 4 64 
and 

n(r..)^^(l + I + i + ^ + ^ + ^ + ...). r^0..^ 0. (230) 
Consequently 

nfr s) 13 1 

-^l-l = l + -r + -r^ + —rs + ---, r ^ 0, s ^ 0, (231) 

where the remaining terms are powers r' s^ of order I + k > 3. 

Similarly, in our study of the birth of a cut we need to estimate K{s) and n(r, s) 
for (r, s) near (1, 1). For K{s) we use the standard formula [27] 

K{s)^-hog{l-s), s-^ 1. (232) 

The appropriate asymptotic expression for n(r, s) is obtained from the inequality 



n r, s] — 



1 + v^(l + v^) VI - V - 
where < r < s < 1 and 



%,.)=tanh-M J^^^ 



< , (233) 



(234) 



To prove (I233p we note that near t = 1 the integrand in (I225p behaves as 
1 1 



X 



vT^{i-rt^)Vi- st^ 72(1 + v^)vr+7^ 

, i , (235) 

v/T^(l - v^t)v/l- v^t 

and that 

dt 20(r, s 



/o v/r^(l- v^t)v/l- v^t v^l - VrVVs-Vr' 
Furthermore, we have 

V2e{r,s) 



(236) 



n(r, s) - 



1 + v^(l + v^) VI - VrVVs - Vr 
' /(r,s,t)-/(r,s,l) 



vr^(i - v^t) vi- 



dt, (237) 



where 
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Now, for all < r, s, t < 1 we have the bound 



so that 



and since 



\dj{r,s,t)\<8 

0</(r,s,t)-/(r,s,l)<8(l-t), 



(239) 
(240) 



< 



< 



(1- v^t)v/l- 
dt 



dt 



(1- v^t)v/l- v^t 
2 



< 



r^/^^y ^/s — ^/r 

7T 



tan 



-1 



,1/4 



s - 



— tan' 



^ ( ^1/4 ^ /T^Ti ^ 



from (123 7p we deduce at once the estimate (12331) . 



(241) 



Appendix C: Asymptotic approximations in the merging and birth of cuts 

In this appendix we consider the two-cut case only and therefore we dispense with the 
superscripts. 

Asymptotic approximations in the merging of two cuts We first prove the asymptotic 
expressions (I172p for (32(T) and /3s(T). Consider the function 

n(r, s) 



C7(/3) = /34-(/34-/33)- 



K{s) 



(242) 



on a neighborhood 

t/ = {/3 G : I A - <6,t = l,...,A} (243) 

of the critical point (3^ = (—2, (3, (3, 2) such that C{(3) is analytic in U. Given (3 E U we 
use (I23ip to approximate C{(3) by its Taylor polynomial at f3^ = {f3i, (32, (^2, f^^) G U: 



C{(3)= (32 + ^(^3-^2) + I ^ 

In particular we have 



C{(3)-(32--{(33-P2) 



<M{(3^-(32)\ (3eU, 



where 



M = - sup 
2 /3e(7 



(245) 



(246) 



Thus, as T — 7- Tf. 



C{(3)^(32 + -{^3-P2) 



(247) 
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Hence, ffT^ and ([MZD lead to 

2 

/^o , (248) 

2 

^ M/33,/3)(/33-/3i)(/33-/34) ' ^^^^^ 
or, equivalently, to 

^(/3i + /34) + ^(3/32 + h) - 2/3) /32 ~ JpTir^y (250) 

^(/3i + /34) + \m + /32) - 2/3) /33 ~ ^^l^- (251) 
If we now set T = Tc — t in fl250p and (125 ip , substitute the approximations 

/3i(r, - t) ~ -2 - /3i(r,)t = -2 + (252) 

/34(T, - t) ~ 2 - P^{T,)t = 2 - (253) 

and look for a solution of the form 

A(Te-t)~/3-^^t'^'+\ -l<ai<0, ^ = 2, 3, (254) 

then fl250D and fl25T|) imply 
1 



\ ^ A fa2+l _ ( ,2a2 + l , ^2 ^3 ga+as+A 

(/3 + 2)2 (/3-2)2; ^ V«2 + l as + l /^^ - 4 



(255) 



(/3 + 2)2 (/3-2)V V«3 + l a2 + l ; /32-4' 

The solution of this system of equations is 

1 1 

02 = CI3 = — , Ao = —A3 = — , (257) 

2' ' 

which completes the proof of (11721) . 

Asymptotic approximations in the birth of a cut In our next discussion all the functions 
involved will be considered as T-dependent functions near Tc. We note that /34 — /Ss — )■ 
and 

l-r~£i^, 1-.^ ^ ^-^-^'^^ (/34-/33), .-r~£i^. (258) 
/3-/3c,2 (^-^(^-M /^-/3c,i 

Thus from (12321) we get 

K{s)r^-^\og{f3,-(3,). (259) 

Let us look for an asymptotic approximation to n(r, s) near T^. Taking into account 
that as T — > 

(/3i,/32,/33,/34)^(^c,iJc,2,A^), (260) 
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/33-/32 



/34-/32 
/34-/3i 



1, 



r 1, 



s-r _ P3-P2 ^ /3-/3c,2 
1 - r /33 - /3i /3 - ' 

a/s-a/t s - r 1 + ^/r P - (5c, 



6'(r, s) — )■ tanh ""^ 



'c,2 



we find that 
Stt 



V2e(r, s 



Therefore it follows from f l233p that 
n(r, s) 



^ 0. 



72 ^(r, s) 



s - -v/r 



Moreover, we have that 



1 - Vr 



1 — r 



s — r 



so that (126 7p can be rewritten as 

tanh-i(,/f^) 

n(r, s) 



Then, using fl258p it follows that 



n(r, s) 



(/3-/3c,i)(/3-/3c,2) 



/94-/33 

which together with (12591) lead to 



tanh ^ 



^/3-/3c,i 



n(r,s) 2^(/3-/3,,i)(/3-/3,,2) 



(/54-/33) log(/34-/33) 



tanh ^ 



c,2 



\ /3 - /3c,i 
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